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ABSTRACT 


A winged midge species, Parochlus steinenii is one of the most abundant species in Antarctica, which is distributed 
over a wide area from the South American continent to the South Shetland Islands in Antarctica. It was dispersed 
into islands in the South Shetland Islands from the South American continent, and it adapted to a variety of environ- 
ments and settled. This species, therefore, is a good model organism to explain the evolutionary process of Antarctic 
terrestrial fauna. Nevertheless, there are few genetic studies on this species, which are necessary for understanding 
the genetic diversity, population structure, etc. Here, we developed and characterized 14 polymorphic microsatellite 
markers. The number of alleles per locus ranged from 2 to 5. The observed and expected heterozygosities were in 
the range of 0.024 to 0.561 and 0.024 to 0.535, respectively. Identifying genetic differences between populations, 
they are suitable markers for researches investigating genetic diversity and population structure of P. steinenii, 
which provide us with clues to dispersion, evolution and ecology of this species. 


Keywords: Parochlus steinenii, Antarctica winged midge, microsatellite marker, next generation sequencing, 


evolution, population genetics 


INTRODUCTION 


Compared with the Arctic that has more than 1,650 species 
of insects, Antarctica has no insects except two dipteran spe- 
cies, winged midge Parochlus steinenii (Gerke) and wing- 
less midge Belgica antarctica Jacobs (Convey and Block, 
1996). Of these, P. steinenii is distributed from the South 
Shetland Islands on the west coast of the Antarctic Peninsula 
to the Tierra del Fuego on the South American continent and 
it is known as one of the most abundant species in Antarctica 
(Brundin, 1970; Convey, 1996). 

Unlike B. antarctica, which lives only in Maritime Ant- 
arctica, P. steinenii was dispersed and settled in several is- 
lands in the South Shetland Islands as well as continent of 
South America (Usher and Edwards, 1984; Edwards and 
Usher, 1985; Convey and Block, 1996). Several researchers 
have suggested that the speciation potential of this species 


which adapted to various environmental conditions may be 
considerable (Sublette and Wirth, 1980; Edwards and Usher, 
1985). This species, therefore, can be a good model to ex- 
plain the distribution and evolution of Antarctic insects. 
But previous studies on P. steinenii have focused mainly on 
morphology and ecology (Wirth and Gressitt, 1967; Edwards 
and Usher, 1985; Shimada et al., 1991; Convey and Block, 
1996; Hahn and Reinhardt, 2006; Rico and Quesada, 2013). 
No suitable genetic markers have been developed for this 
species as well as research to understand the population ge- 
netic structure in order to investigate the evolution of this 
species. 

Population genetics studies are needed to understand the 
adaptation, evolution and dispersion process of organisms. 
One of the most popular genetic markers for population ge- 
netics is a microsatellite due to its high polymorphism, co- 
dominant property, etc. (Tautz, 1989). Information obtained 
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Fig. 1. Map showing the collection sites of Parochlus steinenii in Antarctica. 


using microsatellites is useful for establishing conservation 
measures as well as for investigating the population struc- 
ture, evolution, dispersal and ecological properties (Chistia- 
kov et al., 2006; Selkoe and Toonen, 2006). Here, we firstly 
developed 14 microsatellite markers for P. steinenii. 

Samples of P. steinenii were collected from King George 
Island, West Antarctica (62°14S, 58°47W) in 2015 (Fig. 1, 
site A). Genomic DNAs were extracted using the DNeasy 
Tissue Kit (Qiagen, Valencia, CA, USA). Of these, 41 indi- 
viduals were used in this study. 

Next generation sequencing was performed with the MiSeq 
platform (Illumina, San Diego, CA, USA) and shotgun li- 
braries were prepared and sequenced according to the man- 
ufacturer’s instructions. 

Microsatellite candidates were selected from genomic 
data using the QDD software (Meglécz et al., 2014). A total 
of 504 microsatellite loci were obtained, and for these loci 
56 primers were designed using Primer3 program (http:// 
bioinfo.ut.ee/primer3, Untergasser et al., 2012). PCR am- 
plifications were performed in a total solution of 25 uL con- 
taining 1.0 uL of template DNA, 2.5 uL of 10X buffer, 0.7 
uL of dNTP (10 mM), 1.5 uL of MgCl: (25 mM), 0.5 uL of 
each primer (10 uM), and 0.3 uL of Taq polymerase (Takara, 
Otsu, Japan). Cycling program conditions were follows: ini- 
tial denaturation at 95°C for 5 min, followed by 35 cycles of 
denaturation at 94°C for 60 sec, annealing at 45-60°C for 
60 sec, elongation at 72°C for 90 sec, and final elongation at 
72°C for 7 min. 

M13 tail (FAM: TTT CCC AGT CAC GAC GTT G, PET: 
GCG GAT AAC AAT TTC ACA CAG G, NED: TAA AAC 
GAC GGC CAG TGC) was attached to the 5' end of each 
forward primers, and pig tail(GTT TCT T) was attached to 
the 5’ end of all reverse primers. For genetic diversity anal- 
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ysis, multiplex PCR amplification was performed in 16 uL 
PCR mixture composed of 5 uL of 2x QIAGEN Multiplex 
PCR master mix (Qiagen), 0.08 uL of M13-tailed forward 
primer (10 uM), 0.8 uL of pig tailed reverse primer (10 uM), 
0.3 uL of template DNA, and 0.16 uL of each fluorescence 
primer (FAM/PET/NED) (10 uM). PCR cycle conditions 
were as follows: initial denaturation at 95°C for 15 min, 
90 sec of annealing at a specific temperature (14 cycles at 
63°C, 7 cycles at 58°C, 20 cycles at 55°C), each cycle elon- 
gation at 72°C for 30 sec, and final elongation at 72°C for 
20 min. Genotyping was performed using ABI 3730xl DNA 
analyzer (Applied Biosystems, Foster City, CA, USA) and 
GeneMapper software v.3.7 (Applied Biosystems). 

The number of observed alleles at each locus (A), observ- 
ed heterozygosity (Ho) and expected heterozygosity (He) for 
each locus were calculated using the program ARLEQUIN 
3.5 (Excoffier and Lischer, 2010). Hardy-Weinberg equilib- 
rium (HWE) and linkage disequilibrium (LD) were tested 
using GENEPOP 4.2.1 (Raymond and Rousset, 1995; Rous- 
set, 2008). Program Micro-Checker v2.2.3 was used to 
confirm the existence of null alleles (Van Oosterhout et al., 
2004). 


RESULTS AND DISCUSSION 


A total of 56 primer sets were selected from 504 candidates 
and we successfully developed 14 polymorphic markers 
(Table 1). The number of alleles per locus ranged from 2 to 
5. Ho and Hr varied from 0.024 to 0.561 and 0.024 to 0.535, 
respectively. All loci, except 3 loci, showed significant de- 
viation from HWE after Bonferroni correction (p =0.0035). 
Evidence of null allele was confirmed in 5 loci, which was 
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Table 1. Fourteen microsatellite loci developed from Parochlus steinenii 


Locus Repeat Primer sequence (5'-3') Fluorescent, Size range A Ho He HWE 
motif dye (bp) 

A-10 AAC F: TITCCCAGTCACGACGTTGGCCTTATTTAAAGAATT TAGCAATCG 6FAM 120-123 2 0.079 0.169 0.018 
R: GTTTCTTGTCATGATGGCCTGACCAA 

A-11 AAC F: TITCCCAGTCACGACGTTGTTGTTGGTTAGTGACAACGTCC 6FAM 118-124 2 0.488 0.506 1.000 
R: GTTTCTTAAATTCATAGATGGCTCGAATATC 

A-16 AG F: TTTCCCAGTCACGACGTTGGGTTCCACCGCACTAACACT 6FAM 119-127 2 0.195 0.287 0.069 
R: GTTTCTTGGGCGGAGCCTAAATTTGTA 

B-06 AAC F: TAAAACGACGGCCAGTGCGAGGTGGATTTGTGGCATTC NED 276-279 3 0.075 0.074 1.000 
R: GTTTCTTTCATAGCCGGTGATTTATTCG 

B-11 AAC F: TAAAACGACGGCCAGTGCAACTAACCTGAATTTCGCTAACCA NED 265-268 2 0.244 0.506 0.001? 
R: GTTTCTTGCGCTCAGTTGCCTCAGT 

B-13 AT F: TAAAACGACGGCCAGTGCAAATAAGATGGTGGAGGCGA NED 259-270 5 0.300 0.357 0.372 
R: GTTTCTTGTAAGAAATGTGTATCGGCGG 

B-15 ACCT F: TAAAACGACGGCCAGTGCTGGTGACATTGCTGGAGTTG NED 284-292 3 0.561 0.527 0.603 
R: GTTTCTTCCAACAATATTTGGGCGATT 

B-16 AAG F: TAAAACGACGGCCAGTGCGGCCGTTGTATGACGAAAGT NED 173-183 4 0.128 0.535 0.000° 
R: GTTTCTTTTCATTTCCTTTAATCTTTGAACCA 

C-03 AAC F: GCGGATAACAATTTCACACAGGGGAGAAGTGAGTATTTCGCAGG PET 315-318 2 0.450 0.444 1.000 
R: GTTTCTTCTGTTTGAGTGGTGAAGCTTGT 

C-07 AAG F: GCGGATAACAATTTCACACAGGCAACACCAAATCTTCCTTTGC PET 330-336 4 0.300 0.531 0.001? 
R: GTTTCTTTGCAAATGAATGGCAGAAAG 

C-10 AAT F: GCGGATAACAATTTCACACAGGACCGTTTGAGGATAAAGGAAGA PET 375-378 2 0.289 0.321 0.610 
R: GTTTCTTTTATCCGCTTGCCAATACG 

C-11 AAC F: GCGGATAACAATTTCACACAGGAAATAAATACAGTATCAAGCAGGCA PET 422-423 2 0.024 0.024 1.000 
R: GTTTCTTAGCCCGCCAAGTACTCATT 

C-12 AAG F: GCGGATAACAATTTCACACAGGAGACGCAAATGCTGTGAAAGT PET 429 3 0.024 0.048 0.012 
R: GTTTCTTATCTCACGCCATCACACTGA 

C-13 AAT F: GCGGATAACAATTTCACACAGGGGAAATAGGAGTAGTGCAGTTGG PET 437-439 2 0.171 0.270 0.042 
R: GTTTCTTTCATCTGATCTGGTCAAGGAA 


A, number of allele; Ho, observed heterozygosity; He, expected heterozygosity; HWE, Hardy-Weinberg equilibrium. 


*Significant after Bonferroni correction (p<0.0035). 


presumably due to excess of homozygotes or scoring error. 
The LD of each pair of loci was not significant after Bonfer- 
roni correction. 

The 14 newly developed markers showed low genetic di- 
versity, but it was probably because the samples used for the 
analyses were collected from an isolated area. Although not 
presented in this paper, we performed Bayesian clustering of 
this population with the population of Deception Island (Fig. 
1, site B) in the South Shetland Islands using the STRUC- 
TURE program (Pritchard et al., 2000). As a result, we con- 
firmed that the two populations had a genetically distinct 
structure. These results confirmed that the newly developed 
markers show discernment to investigate the genetic struc- 
ture of populations. These new markers could be useful for 
establishing a conservation management plan of indigenous 
species in Maritime Antarctica. 
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